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Abstract 



> 

X. 

Inclusion of the QED higher order radiative corrections in the two- 
photon process, e + e~ — > e + e~fj> + fi~, is examined by means of the 
structure function and the parton shower. Results are compared with 
the exact 0(a) calculations and give a good agreement. These two 
methods should be universally applicable to any other non-annihilation 
processes like the single- W productions in the e + e~ collisions. In this 
case, however, the energy scale for the evolution by the renormalization- 
group equation should be chosen properly depending on the dominant 
diagrams for the given process. A method to find the most suitable 
energy scale is proposed. 



1 Introduction 



A precise prediction of the cross sections for high-energy e + e~ scattering fre- 
quently requires an estimation of the radiative corrections beyond the lowest 
order calculations. Among various corrections from the electro-weak inter- 
actions it is known that the QED radiative corrections to the initial-state 
particles gives the greatest contribution in general. Sometimes all order sum- 
mation is necessary to give the needed precision. In the leading-logarithmic 
approximation the higher-order summation of the QED corrections can be 
done easily thanks to the factorization theorem. The theorem given by ref . [|IJ 
guarantees that the final state corrections result only a small contribution to 
the total cross section. 

For the e + e~ annihilation processes we have tools. They are universally 
applicable because only the initial-state QED corrections are involved; the 
structure function || and the parton shower H, |[] methods widely used in high 
energy physics today. The recent experiments, however, require the higher 
order corrections for the multi-particle final states, such as the four-fermion 
productions at LEP2. Even though the exact calculations of the higher order 
corrections are very difficult or impossible, still it is possible to include the 
biggest QED corrections by making use of those tools as long as the process 
is dominated by the annihilation. 

On the other hand for the non- annihilation processes it has been not well 
investigated how to apply these tools. Only a few examples are the Bhabha 
scattering to which the structure function has been used in|| and the parton 
shower in||. In the present work it will be shown that the structure function 
and the parton shower can be also the universal tools for any non-annihilation 
process once the evolution energy is settled. 

First we apply and examine both methods to the two-photon process, 
e + e~ — > e + e~fi + fi~ in the next section. The obtained total and differential 
cross sections are compared with those given by the 0(a) corrections |J in 
section 3. These methods must be universal to any process, again thanks to 
the factorization theorem. The energy scale, however, with which the radia- 
tive correction is evolved should be carefully chosen for the non-annihilation 
processes. An unique and definite way to find this energy scale is explained 
in section 4. 
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2 Calculation method 



2.1 Structure Function Method 



The structure function(SF) for the initial state radiative correction(ISR) is 
well established for the case of the e + e~ annihilation processes. The observed 
cross section corrected by ISR can be expressed by using SF as 

vtotai(s) = J dx 1 J dx 2 D e -(x 1 ,s)D e +(x2,s)ao(x 1 x 2 s), (1) 

where D e ±(x,s) is the electron(positron) structure function with xi }2 being 
the energy fractions of e^. The SF up to the 0(a 2 ) is given byj7| 
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In deriving this formula, which is given by solving the Altarelli-Parisi equa- 
tion in the LL approximation^, we have used an ad hoc trick to get more ac- 
curacy. The factor f3 = {2a / ti)1h{s / m 2 e ) is replaced by (3 = (2a/7r)(ln(s/mg) — 
1) to match with the perturbative calculations. 

Let us apply the SF method to the two-photon process, 



e (p_) +e + (p + ) 



e (g_) + e+(g + )+/i-(£;_) +/!+(£;+). 



(4) 



For the forward scattering of e^, the multi-peripheral diagrams shown in 
Fig.l give the dominant contribution to the total cross section. Thus only 
the multi-peripheral diagrams are taken into account in this work. In this 
case the corrected cross section is given by 

o- t otai(s,t±) = J dxj- J dx F - J dx I+ J dx F+ D e -(x I -,Q 2 _)D e -(x F ^,Q 2 _) 
D e +(x I+ , Q 2 + )D e +(x F+ , Q 2 + )a (s, t±). (5) 

Here Q\ is the energy scale to be fixed, with which SF should be driven. 
Since these functions are common to the initial and the final radiations from 



2 



the e s in the leading-logarithmic(LL) approximation, we shall drop the 
subscript from the SF hereafter. After(before) the photon radiation the 
initial(fmal) momenta p± (q±) become p± (q±) in the following ways 

P- = Xi-P-, q- = ~^—q-, (6) 

x F ^ 

P+ = xi+P+, Q+ = —<1+, (7) 

x F+ 

respectively. Then the CM energy squared(s = (p_ + p+) 2 ) and the momen- 
tum transfer squared(t± = (p± — q±) 2 ) are scaled as follows, 

Xj± 

s = xi-Xi + s, t± = 1±. (8) 

x F± 

Note that SF behaves like 6(1 — x) when (3 — > 0, that is, 

ix^-{l-x)i- l f{x) = f(l)+ f 1 dx\(l-x)l -l]f(x) 
2 jo L J 

= /(l) + t f 1 dxln(l - x)f'(x) + 0((3 2 ), (9) 



2 Jo 

where f(x) is an arbitrary smooth function. 

The choice of the energy scale in SF is not a trivial matter. As pointed out 
in Ref . [[llj it is natural to use —t± instead of s for the two-photon processes. 



The justification of this choice is given by comparing them with the pertur- 
bative calculations. This could be done in the region where the soft-photon 
approximation of the total cross section is valid. Since the corrections for the 
e~ and e + sides must be symmetric, it is enough to consider only those from 
the e~ side. The total correction will be obtained by doubling them. Then 
Eq.(||) can be simplified as 

cr(s) = J dxj J dxFD(xi,Q 2 )D(xF,Q 2 )<Jo(xis). (10) 

The integrations are performed in the region 1 — xi t p 1. Since the Born 
cross section of this process(<T (s)) is a smooth function of s, the cross section 
in the soft photon approximation is written as 

Vsojt = cro(s) f dxi [ dx F D(x I ,Q 2 )D(x F ,Q 2 ) 
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where k c is the maximum energy of the sum of the initial and the final 
photon energies and E = p° « q°. The function H is called radiator and can 
be obtained easily from SF as 
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The integration of the function H in the small k c / E region gives 

3 
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where L = ln(Q 2 /ml), and / = h\.(E/k c ). Then the cross section in the soft 
photon approximation up to the 0(a) is obtained, 
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This expression is compared with the perturbative calculation given by Berends, 
Daverveldt and Kleiss0(BDK hereafter). In the BDK program the multi- 
peripheral diagrams and their 0(a) corrections of the self-energy correction 
to e^, the vertex correction for e ± -e ± -7 vertex, the soft- and hard-photon 
emissions and the vacuum polarization of the virtual photons are calculated. 
The corrections from the photon bridged between different charged lines are 
not given because the contributions from the box diagrams with photon ex- 
change between e + and e~ is known to be small [|lOfl . The LL approximation 
of the virtual correction factors (vertex + soft photon) is 



2ReF 1 + 5 S 
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where L t = ln(—t/m 2 ) and t = (p_ — q-) 2 . By comparing Eqs. (TT3J) and (|T6|) 
one concludes that the energy scale of SF should be Q 2 = t 

This comparison shows that once the proper energy scale is found the 
SF can reproduce the correct evolution of the soft-photon emission along the 
electron line up to the 0(a). However, there remains some mismatch with 
the constant term. This can be compensated by multiplying an overall factor 
to the SF. This factor is usually called K-factor. When the SF is evolved 
with —t this factor is found to be 1 — a/2ir. Then the total K-factor for both 
e + and e~ must be 1 — a/n. 

It should be noted that the assumption —t/m 2 3> 1 does not hold for the 
forward scattering of the two-photon process. The LL-terms in the SF are 
no more leading if this happens. To find the region where the approximation 
is valid the LL terms of the form factor Fx is compared with the exact form 
given by BDK in Fig. 2. While they agree well at high Q 2 region as expected, 
they show a large deviation in the region Q 2 /m 2 e < 10. In order to make 
SF well-defined, the energy evolution must be truncated at some point, say 
L = 1. 

2.2 Parton Shower Method 

Instead of the analytic formula of the structure function, a Monte Cairo 
method based on the parton shower algorithm in QED (QEDPS) can be used 
to solve the Altarelli-Parisi equation in the LL approximation^]. The de- 
tailed algorithm of the QEDPS is found in Ref.[Q for the e + e~ annihilation 
processes and in Ref. || for the Bhabha process. For the two-photon process 
the energy scale for the parton shower evolution is also —t. The truncation 
at L — 1 is again imposed. One difference between SF and QEDPS is that 
the ad hoc replacement L — > L — 1, which was realized by hand for SF, can- 
not be done for QEDPS. This causes a deviation of the K-factor from the SF 
method. The K-factor for QEDPS is given by 1 — 4a /ir. Another significant 
difference between these two is that QEDPS can treat the transverse momen- 
tum of emitted photons correctly by imposing the exact kinematics at the 
e — > cy splitting. It does not affect the total cross sections so much when the 
final e ± have no cut. However, the finite recoiling of the final may result 
a large effect on the tagged cross sections as shown below. 

In return to the the exact kinematics at the e — > cy splitting, e ± are no 
more on-shell after photon emission. On the other hand the matrix element 
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of the hard scattering process must be calculated with the on-shell external 
particles. A trick to map the off-shell four-momenta of the initial to those 
at on-shell is needed. Following method is used in the calculations. 

1. s = (p_ + p + ) 2 is calculated, where p± are the four- momenta of the 
initial e after the photon emission by QEDPS. s can be positive even 
for the off-shell e ± . (When s is negative, that event is discarded.) 

2. All four- momenta are generated in the rest-frame of the initial e after 
the photon emission. Four-momenta of the initial e in this frame are 
p±, where p± = ml (on-shell) and s = (p_ + p+) 2 

3. All four-momenta are rotated and boosted to match the three-momenta 
of p± with those of p±. 

This method respects the direction of the final e ± rather than the CM energy 
of the collision. The total energy does not conserve because of the virtuality 
of the initial e ± . 



3 Numerical Calculations 



3.1 No-cut Case 



First the total cross section of e + e — > e + e without any experimen- 

tal cut is considered. The exact matrix element is generated by the GRACE 



system 1 12[|. Only the multi-peripheral diagrams are generated for the test 



of the approximation methods. The phase-space integration of the matrix 



element squared is carried out numerically by BASES |L3| using an adaptive 
Monte Cairo method. It is trivially confirmed that the Born cross sections 
agree with BDK results within the statistical error of the numerical integration. 

The convolution of SF or QEDPS with the Born cross section is not a 
straightforward task because of the complicated four-body kinematics. For 
a simple two-body process such as the Bhabha scattering the momentum 
transfer (— t) can be taken as one of the integration variables and the Eq.(^) 
will be easily performed then. It is, however, practically impossible to use 
this technique for the case of the four-body kinematics. Hence we have to 
admit the following approximation: 



6 



1. Eight random numbers obtained from BASES correspond to the eight 
independent variables of the phase-space integration. 

2. All the kinematical variables are determined with no radiative effect. 
The — £±'s are fixed also. 

3. The radiative correction factors of Eq.(5) are determined from — £±'s, 
x/_ and xj + . 

4. A new CM energy is obtained by the fixed xi- and xi + . 

5. All the kinematical variables are re-calculated from the same set of the 
random numbers and the new CM energy. 

In this method the evolution energy scale of the radiator is not exactly the 
same after the re-calculation with the new CM-energy. However, this differ- 
ence must be beyond the LL order. 

The total cross sections without any experimental cuts with SF and QEDPS 
are summarized in Table |l|. The BDK program also includes the correction 
from the vacuum polarization of the photon propagator. In order to compare 
the SF and QEDPS results with BDK, this correction is removed. After including 
the K-f actor the results of our two methods are in good agreement with the 
0(a) calculation(without the vacuum polarization). 



E CM (GeV) 


0o {nb) 


a sf (nb) 


oqedps (nb) 


(Jbdk {nb) 


20 


97.0(1) 


96.3(1) 


96.0(2) 


96.0(1) 


40 


137.5(3) 


136.3(1) 


135.9(3) 


135.9(1) 


100 


202.8(4) 


201.0(2) 


200.5(4) 


200.5(2) 


200 


262.0(6) 


259.8(3) 


259.1(6) 


258.8(2) 



Table 1: Total cross sections without the vacuum correction. Co is the tree cross section. 
Last column shows the results by the BDK program with the vertex correction, the soft- 
photon correction and the hard-photon emission. The vacuum polarization is removed from 
BDK. The number in parenthesis shows the statistical error of the numerical integration on 
the last digit. 

The effect of the vacuum polarization can be included if one uses the 
running QED coupling in the SF method. The ffj coupling evolved by the 
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renormalization group equation is given by 

^/7H ± )=^/7(0)(l-^E^log^e(-t ± -m^ , (17) 

where a = 1/137.036 is the QED coupling at the zero momentum transfer, 
Ci the color factor, the electric charge in unit of the e + charge, the mass 
of the i-th fermion. The index i runs over all massive fermions. Only those 
fermions whose mass is greater than —t± are taken into account through the 
step function. The quark masses are chosen so as to match with the vacuum 
polarization in the BDK program. The results are shown in Table ^. The 
deviation from the BDK program is typically around 0.5%. 



Ecu (GeV) 


cr sf {nb) 


oqedps {nb) 


(Jbdk {nb) 


20 


97.5 (1) 


96.9(1) 


97.1(1) 


40 


137.8 (1) 


137.1(2) 


137.3(1) 


100 


203.3 (1) 


202.2(3) 


202.4(2) 


200 


262.6 (1) 


261.2(3) 


261.2(2) 



Table 2: Total cross sections with the vacuum polarization. The number in parenthesis 
shows the statistical error of the numerical integration on the last digit. 

The differential cross sections with respect to the e~ energy and angle, 
the CM energy of the final four-fermions and the invariant mass of the fi^- 
pair at the CM energy of 200 GeV are shown in Fig. 3. The two programs 
give consistent distributions. 

In order to check the recoiling effect of the final e~ due to the photon 
emission by QEDPS, the e~ polar angle is compared in Fig.4 between SF and 
QEDPS. The cross sections with the e~ angle between 10° and 20° are found 
to be 6.00(5.80) pb for QEDPS(SF) at the CM energy of 200 GeV. It is not 
surprising that the agreement becomes worse than the total cross sections, 
because SF includes no recoiling at all. 

If a wrong energy scale of s = (p_ + p + ) 2 is used instead of t± as the 
energy evolution scale in the ISR tool, one may get over-estimation of the 
ISR effect. At the CM energy of 200 GeV, SF with the energy scale of s gives 
the total cross section of 257.8nb instead of 262. 6nb with the correct energy 
scale. 
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3.2 Single-tagging Case 

The same comparison is done for the e 
applied are: 
For the e~, 

• io° < e e - < 170°, 



-tagging case. The experimental cuts 



• E e - > 1 GeV. 
For the /i ± , 

• 10° < M ± < 170°, 

• E IX ±>1 GeV, 

• M m > 1 GeV. 

The total cross sections with the above cuts at the CM-energy of 200 GeV 
are calculated to be 1.169 ± 0.004pb(l. 13 ± O.Olpb) by GRACE with QEDPS(by 
BDK). The vacuum polarization, i.e. the running a, is included. This small 
discrepancy may come from the finite recoiling of the final e~ by the soft 
photon emissions in QEDPS. The differential cross sections are also compared 
in Fig. 4. The results of GRACE with QEDPS are in good agreement with BDK. 



4 Energy Scale Determination 

The factorization theorem for the QED radiative corrections in the LL ap- 
proximation is valid independent of the structure of the matrix element of 
the kernel process. Hence SF and QEDPS must be applicable to any e + e~ scat- 
tering processes. However, the choice of the energy scale in SF and QEDPS 
is not a trivial issue. For a simple process like the two-photon process with 
only the multi-peripheral diagrams considered so far, the evolution energy 
scale could be determined by making use of the exact perturbative calcula- 
tions. However, this is not always possible for more complicated processes. 
Hence a way to find a suitable energy scale without knowing the exact loop 
calculations should be established somehow. 

First let us look at the general consequence of the soft photon approxima- 
tion. The soft photon cross section(including both the real and the virtual 
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photon effects) is given by the Born cross section multiplied by some correc- 
tion factor in the LL order as [14] 



da, 



oft{ 



dtt 



x 



Pi 




(18) 



(19) 



where m/s^'s) are the mass(momentum) of j-th charged particle, k c the 
maximum energy of the soft photon (boundary between soft- and hard- 
photons), E the beam energy, and ej the electric charge in unit of the e + 
charge. The factor rjj is —1 for the initial particles and +1 for the final par- 
ticles. The indices run over all the charged particles in the initial and 
final states. 

For the process (f|) one can see that the soft-photon factor in Eq.flTSP 
with a (p_ ■ g_)-term reproduces Eq. ([T6|) in the LL approximation. This 
implies that one is able to read off the possible evolution energy scale in SF 
from Eq. ([Tj|) without explicit loop calculations. However, one may have a 
question why the energy scale s = (p_ + P+) 2 does not appear in the soft- 
photon correction even they are included in Eq.([H|). When we applied SF to 
the two-photon process in the previous section, we have ignored those terms 
which come from the photon bridged between different charged lines. This 
is because the contributions from the box diagrams with photon exchange 
between e + and e~ is known to be small fT0|] . Fortunately the infrared part of 
the loop correction is already included in Eq. fll8|) and no need to know the full 
form of the loop diagram. For the two-photon processes if one looks at two 
terms with, for example (p_ ■ p+)- and with (g_ ■ p + )-terms, the momentum 
of e~ is almost the same before and after the scattering(p_ « q_). Only 
the difference appears in rjjTjk = +1 for a (p_p + )-term and T)jT]k = —1 for a 
(g_p + )-term. Then these terms compensate each other after summing them 
up for the forward scattering which is the dominant kinematical region of 
this process. This is why the energy scale s = (p_ +p+) 2 does not appear in 



the soft-photon correction despite the fact that it exists in Eq. (|18"D . 

When some experimental cuts are imposed, for example the final e~ is 
tagged in a large angle, this cancellation is not perfect but partial and the 
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energy scale s must appear in the soft-photon correction. In this case the 
annihilation type diagrams will also contribute to the matrix elements. Then 
the usual SF and QEDPS for the annihilation processes are justified to be 
used for the ISR with the energy scale s. One can check which energy scale 
is dominant under the given experimental cuts by numerically integrating 
the soft-photon cross section given by Eq. (|T8|) over the allowed kinematical 
region. Thus in order to determine the energy scale it is sufficient to know 
the infrared behavior of the radiative process using the soft-photon factor. 
In some region of the phase-space two or more energy scales may be involved 
in the soft-photon cross section with comparable amount of contribution. In 
this region a simple SF and QEDPS are not applicable. 

5 Conclusions 

Two practical tools to incorporate with the QED radiative corrections are 
developed for non-annihilation processes by means of the structure function 
and the parton shower. These programs are applied to the two-photon pro- 
cess, e + e~ — ► e + e~/i + yU~. The results are compared with the perturbative 
calculation of the 0(a) and show a good agreement. 

These tools should be applicable to any non-annihilation process uni- 
versally. It is demonstrated that the energy scale for the evolution, which 
depends on the dominant diagrams in the interested kinematical region, can 
be determined with the help of the well known formula of the soft photon 
factor. As an example we have tried a real W production in e + e~ anni- 
hilation. The application of these tools to more complicated processes like 
four-fermion final state including a single- W production are left to the future 
publications. 
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Figure 1: Multi peripheral diagrams. 




Figure 2: The exact form- factor (Rei^) and that with LL approximation. 
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Figure 3: The differential cross sections without any experimental cuts. The solid his- 
tograms show the GRACE-with-QEDPS results, the stars the GRACE-with-SF, and the circles 
the BDK results. 
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Fi gure 4: The e angle distribution with respect to the beam axis without any exper- 
imental cuts. The solid histograms show the GRACE-with-QEDPS results, the stars the 
GRACE-with-SF, and the circles the BDK results. 
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Figure 5: The differential cross sections with the e -tagging and visible muon-pair. The 
solid histograms (circles) show GRACE with QEDPS (BDK) results. 
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